Weyl Fermions and the Anomalous Hall Effect in Metallic Ferromagnets 
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We reconsider the problem of the anomalous Hall effect in ferromagnetic SrRuOs, incorporating 
insights from the recently developed theory of Weyl semimetals. We demonstrate that SrRuOs 
possesses a large number of Weyl nodes, separated in momentum space, in its bandstructure. While 
the nodes normally do not coincide with the Fermi energy, unless the material is doped, we show 
that even the nodes inside the Fermi sea have a significant effect on the physical properties of 
the material. In particular, we show that the common belief that (non-quantized part of) the 
intrinsic anomalous Hall conductivity of a ferromagnetic metal is entirely a Fermi surface property, 
is incorrect: there generally exists a contribution to the anomalous Hall conductivity that arises 
from topological Fermi-arc surface states, associated with the Weyl nodes, which is of the same 
order of magnitude as the Fermi surface contribution. 



Understanding topological properties of the electronic 
structure of materials and their experimentally observ- 
able consequences has become one of the most impor- 
tant themes of the modern condensed matter physics, 
gradually replacing the traditional focus on the proper- 
ties, determined by symmetry. This shift has become 
particularly obvious in recent years, after the discovery 
of topological insulators pTj. An interesting new devel- 
opment in this field is the realization that to possess 
topologically nontrivial properties a material does not 
have to be an insulator, which significantly expands the 
range of potential realizations of such materials [2H7J • In 
particular, it is now understood that nondegenerate ac- 
cidental band crossings, named Weyl nodes in Ref. [3], 
which occur generically in any three-dimensional mate- 
rial with a broken time reversal or inversion symmetry, 
have topologically nontrivial properties. These crossings 
act as "magnetic monopoles" in momentum space and 
the corresponding quantized topological charge gives the 
band-crossing points stability to perturbations. 

An alternative view of the Weyl nodes, which ap- 
plies specifically to ferromagnetic materials, in which 
the nodes appear due to broken time reversal symmetry, 
is based on regarding the three-dimensional (3D) band- 
structure as a set of two-dimensional (2D) bandstruc- 
tures, parametrized by the crystal momentum compo- 
nent k z , where the z-axis is taken to be along the mag- 
netization direction. Weyl nodes in this picture corre- 
spond to gap-closing quantum phase transitions, at which 
Chern numbers of the 2D bands change pairwise by plus 
and minus the topological charge of the corresponding 
node. In other words, a Weyl node can be viewed as a 
quantum Hall transition point in momentum space. One 
interesting consequence, that immediately follows from 
this, is the existence of chiral Fermi-arc surface states, 
corresponding to pairs of bands with nonzero Chern num- 
bers, where the Fermi arc connects projections of the 
Weyl nodes with opposite topological charge on the sam- 



ple surface Brillouin zone (BZ) [31 [5]. 

As we show in this paper, this viewpoint is particu- 
larly useful in discussing the role of the Weyl nodes in 
intrinsic anomalous Hall effect (AHE) in metallic ferro- 
magnets. AHE in ferromagnets is an old problem [5], 
the interest in which has been revived recently with the 
realization that topological properties of the electronic 
structure likely play a very important role in it [9lH6j. 
In particular, here we demonstrate that Weyl nodes, and 
the associated topological Fermi-arc surface states con- 
tribute significantly to the intrinsic anomalous Hall con- 
ductivity of ferromagnetic metals, which implies that the 
anomalous Hall conductivity in ferromagnets can not be 
viewed as a Fermi surface property. 

While our results likely apply to most metallic ferro- 
magnets, we choose SrRuOs as the specific material we 
focus on. The main reason for this choice is a relative sim- 
plicity of its relevant electronic structure near the Fermi 
energy, which to a good approximation consists of six 
bands, corresponding to the three t 2g Ru d-orbitals [17] , 
which are 2/3-fillcd in the undoped material. We also as- 
sume undistortcd cubic perovskite crystal structure for 
SrRuC>3, which should be a good approximation, espe- 
cially since we are interested in topological properties of 
the bandstructure, insensitive to minor variations of the 
parameters. 

We start from a 6-orbital tight-binding model for 
SrRuC>3, which can be written down based on symmetry 
considerations |18) . The momentum space Hamiltonian 
has the following form: 
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where summation over repeated orbital and spin indices 
is implicit. The first term in Eq. (nl) corresponds to spin- 
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split unhybridized t2g-orbital band dispersions: 
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where t\ and £2 are the in-plane and out-of-plane or- 
bital hopping matrix elements, m is the rigid momentum- 
independent exchange spin splitting, which is nonzero in 
the ferromagnetic phase, and crystal momentum is mea- 
sured in units of the cubic structure lattice constant. Sec- 
ond term corresponds to inter-orbital spin-independent 
hopping, which is induced by hybridization between the 
oxygen p-orbitals and ruthenium d-orbitals: 



f£ b = -4/sin(fc Q )sin(fc b ). 
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Finally, the third term corresponds to spin-orbit (SO) in- 
teractions, projected onto the tig orbital manifold, where 
A denotes the strength of the SO interactions, e abc is the 
fully antisymmetric tensor and r 1 ' 2,3 is the triplet of Pauli 
matrices. We should point out here that we will not at- 
tempt any quantitative comparison with experiments in 
this paper and treat parameters in Eq. ([I]) with some 
freedom. What we are after in this work is qualitative 
features of the bandstructure and their role in the AHE, 
and we expect our results to apply quite generally to 
metallic ferromagnets, not just to SrRu03. 

In the absence of the SO interactions, the eigenstates 
of H exhibit multiple band crossings, which arise due to 
the spin splitting. One might think that once the SO 
interactions are turned on, these band crossings will be 
eliminated and replaced by gaps of magnitude, propor- 
tional to the SO interaction strength A |TB] for small A. 
Recent work on Weyl semimetals, however, has demon- 
strated that this is not the case: nondegenerate band 
crossings in three dimensions possess topological stabil- 
ity and can not be eliminated by small perturbations. 
This is illustrated in Fig. [T] We find that a large number 
of band crossings in SrRu03 survive even in the presence 
of the SO interactions and have a very significant quali- 
tative effect on the anomalous Hall conductivity, as will 
be shown below. 

As is well-known, intrinsic anomalous Hall conductiv- 
ity can be calculated as an integral of the Berry curvature 
over all the occupied states |16j : 
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where np(e) is the Fermi-Dirac distribution function, 
e n k is the band dispersion, and l~2 ri k is the Berry cur- 
vature vector, which is a curl of the Berry connection 
A„k = «(nk|Vfe|nk), and may be thought of as an ana- 
log of magnetic field in momentum space. One way to 



1 




^3 -1 




-— — ^7 





(a) 
EAi 



z 





















(b) 



FIG. 1. (Color online), (a) Bandstructure of SrRuC>3 along 
the k x = k y = line with \/U = 0,t 2 /ti = -0.2, f/h = 
0.2, and m/ti = 1. Nonzero m induces crossings between 
several of the bands, (b) Same bandstructure but with A/ti = 
0.4. While some of the crossings are eliminated, most clearly 
survive. 



understand the important role, played by the Weyl nodes 
in the intrinsic AHE, is to realize that these nodes act 
as magnetic-monopole-like sources of the Berry curva- 
ture field, which is divergenceless in the absence of such 
sources. For our purposes, however, it will be more useful 
to adopt a different viewpoint. Namely, since the magne- 
tization m along the z-axis introduces a preferred direc- 
tion and reduces the cubic symmetry down to tetragonal, 
we can separate the 3D momentum space integration in 
Eq. Q into a ID integral over k z and a 2D integral over 
kx = (k x , ky): 
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where 
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We may then regard a X y(k z ) as the Hall conductivity 
of a set of 2D systems, parametrized by k z . This repre- 
sentation of o X y turns out to be extremely useful, as we 
show below. 
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FIG. 2. (Color online). Band dispersion for fe/ti = 
-0.2, f/ti = 0.2, A/ti = 0.4,m/ti = 1, plotted along high- 
symmetry directions using tetragonal lattice system notation. 
There are multiple band crossings between different bands. 
Some of the crossings are in fact away from high-symmetry 
directions and can not be seen in this plot. 

As was first pointed out by Haldane [H], we can use 
Stokes theorem and rewrite a 2 ® in the following way: 

C(y = TrE/ dk ' A * (7) 

n n J 

where the integral is along the ID Fermi surface of a 
2D system, corresponding to a given k z . This suggests 
that the intrinsic anomalous Hall conductivity may be 
thought of as a Fermi surface property, like all other 
transport properties of metals [14] . However, there is 
a subtlety in this argument. The Stokes theorem applies 
only when the corresponding band has a zero Chern num- 
ber [Tp , 20 . Chern number is an obstruction to the Stokes 
theorem [20] and the correct form of Eq. ^ is: 
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where 

c n {k z ) = ± J d 2 k ± n* lk jk z ), (9) 

is the Chern number of the two-dimensional band n at 
momentum k z and the sum over n in the second term in 
Eq. p| is restricted to completely filled 2D bands only, 
which is indicated by the prime. To make Eq. ^ well- 
defined we may regard the first BZ as an open square 
rather than a torus, so that when a band is completely 
filled, the line integral in the first term becomes an inte- 
gral over the BZ boundary and gives the corresponding 
Chern number. Gauge ambiguity in the first term may 
be eliminated by comparing Eq. ^ with the explicitly 
gauge- invariant result, obtained from Eq. (§. 

If Weyl nodes were not present in the bandstructure, 
all the Chern numbers C n (k z ) would not in fact depend 
on k z , and it is normally assumed that they are zero, 
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FIG. 3. (Color online), (a) Chern number for the second low- 
est band as a function of k z for the same parameter values as 
in Figs. [T] [2] The Chern number experiences multiple integer- 
valued jumps as a function of k z in the first BZ. The jumps 
are the evidence of the existence of Weyl nodes, (b) Locations 
of all Weyl nodes in the bandstructure of SrRuOs. Only one 
quarter of the first Brillouin zone is shown. The remaining 
points are obtained by successive 7r/2-rotations around the 
z-axis. Dots are linear (charge one) Weyl nodes, while crosses 
denote quadratic (charge two) Weyl nodes. The solid (con- 
necting pairs of nodes with opposite topological charge) and 
dashed lines are guides to the eye. 



since a very large spin splitting would be needed to cre- 
ate a band with a constant nonzero C n (k z ), implying 
a contribution to the total Hall conductivity from this 
band, quantized in units of e 2 G/2irh, where G — 2ir is 
the smallest reciprocal lattice vector. Thus one may con- 
clude that the anomalous Hall conductivity is entirely a 
Fermi surface property |14j . This, however, is incorrect, 
since C n (k z ) are not independent of k z in the presence 
of Weyl nodes, and they arc in general present in any 
ferromagnet, even when the magnetization is small. In- 
deed, as discussed above, a Weyl node may be thought 
of as a gap-closing transition point for the 2D bandstruc- 
ture, parametrized by k z . At the transition point, the 
Chern numbers of the two bands, that touch at the Weyl 
node, change by ± the topological charge of the node [3]- 
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6 . The corresponding contribution to the total 3D Hall 
conductivity is then not quantized in units of G, but is 
still not a Fermi surface property: it instead corresponds 
to the Fermi-arc surface states, associated with pairs of 
Weyl nodes with opposite topological charge. 

Let us now see how the above picture is realized 
in SrRuOa. The 6-band tight-binding Hamiltonian of 
SrRu03 Eq. is easily diagonalized numerically. The 
bandstructure, plotted along high-symmetry directions 
in the first BZ in the standard way, is shown in Fig. [2] 
Some of the Weyl nodes are in fact clearly visible in 
Fig. [2j In general, however, finding all the nodes, along 
with their topological charges, is quite a difficult task. 
What makes it significantly easier is the relation between 
the Weyl nodes and the change of the Chern number as 
a function of k z , discussed above. To evaluate the Chern 
numbers, we use a discrete lattice version of Eq. (9) [25, 
corresponding to the discrete values of the crystal mo- 
mentum in a finite sample with periodic boundary con- 
ditions. Namely we calculate the lattice Berry connection 
fields as: 
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where /i are the nearest-neighbor vectors of the square 
momentum-space lattice with a fixed k z . The Chern 
number is then evaluated as: 



C n (kz) 
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An example of the resulting plot of C n (k z ) is shown in 
Fig. [3] Every discrete jump of the Chern number by an 
integer is due to the presence of Weyl nodes in the plane 
in 3D BZ, corresponding to the given value of k z . The 
magnitude of the jump is equal to the total topological 
charge of Weyl nodes in a given plane. The informa- 
tion obtained from these plots, plus symmetry consider- 
ations, confirmed by an explicit examination of the band 
dispersions, allow us to identify all the Weyl nodes in 
the bandstructure of SrRuC>3, as shown in Fig. [3j We 
find a total of 22 pairs of nodes with opposite topolog- 
ical charges in the first BZ. Two of those pairs, sepa- 
rated along the (k x ,k y ) = (0,0) and (k x ,k y ) = (ir,ir) 
lines, correspond to double- Weyl nodes with topological 
charges ±2. These are stabilized by the cubic symmetry 
of the undistorted perovskite structure and will be split 
each into two pairs of ordinary Weyl nodes with single 
unit of topological charge when the orthorhombic distor- 
tion, typically present in real SrRu.03 material, is taken 
into account [22] ■ The rest of the nodes have topolog- 
ical charge ±1. Of those, four pairs are located on the 
same lines as the double- Weyl nodes, four pairs are on 
the (k x ,k y ) = (0, 7r), (it, 0) lines, eight pairs are located 
in the k x — ±k y planes, and four pairs are located in the 
k x = 7r and k y = tt planes. The nodes are generally not 
at the same energy, unless required by symmetry. 




00 



FIG. 4. (Color online). Total anomalous Hall conductivity 
T xy (solid line) and the edge-state contribution to the con- 



ductivity a^y 9e (dashed line) as a function of Fermi energy ef 
for m/t 1 = 1 (a), and magnetization m for e^/ii = (b) 
is in units of e 2 /h (the lattice constant is set to unity). 
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We can now explicitly identify two distinct contri- 
butions to the total anomalous Hall conductivity of 
SrRuC>3: the contribution arising from the Femi-arc sur- 
face states and the contribution, associated with the bulk 
Fermi surface. We evaluate the surface state part of the 
Hall conductivity by summing the Chern-number contri- 
butions of all completely filled 2D bands at every k z [5U] : 
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(12) 



The bulk Fermi surface contribution is then the difference 
between the total Hall conductivity, given by Eq. Q, 
and the surface state contribution, Eq. (12). The total 



anomalous Hall conductivity is evaluated by summing 



the lattice Berry curvature, Eq. ( 11 ), over all states below 



the Fermi energy. The results are shown in Fig. [4j It is 
seen that a xy is a highly nonmonotonic function of both 
the Fermi energy (which can be varied to some degree 
by doping) and the magnetization. The peaks and dips 
in the dependence of a xy onej? correspond to the Fermi 
level passing through the Weyl nodes. It can also be seen 
that the Fermi-arc surface state contribution is always of 
the same order as the Fermi surface contribution and can 
not be regarded as an insignificant correction. 

Eq. ( 12 ) appears to suggest that all Weyl nodes be- 
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low the Fermi energy contribute to the anomalous Hall 
conductivity. This, however, is not correct. As discussed 
above, Weyl nodes may be thought of as points of quan- 
tum Hall transitions in momentum space, at which the 
corresponding 2D bands experience equal in magnitude 
but opposite in sign change of the Chern number. This 
means, in particular, that when both 2D bands in the 
pair, joined by a pair of Weyl nodes, are filled, the total 
contribution of this pair of nodes to the anomalous Hall 
conductivity is zero. 

In conclusion, we have identified the role, played by 
nondegenerate band-touching nodes (Weyl nodes) in the 
intrinsic AHE in ferromagnetic metals. We have shown 
that, in general, even the non-quantized part of the 
anomalous Hall conductivity is not a Fermi liquid prop- 
erty, in the sense that a significant part of it arises from 
the Femi-arc surface states, associated with pairs of Weyl 
nodes of opposite chirality, and not with bulk states on 
the Fermi surface. It is clear from our analysis that any 
generic model of intrinsic AHE in ferromagnetic met- 
als must incorporate Weyl nodes and the currently used 
models are lacking in this regard. 

While the Fermi surface and Fermi-arc edge state con- 
tributions to the anomalous Hall conductivity are clearly 
separated in theory, experimentally one of course mea- 
sures the total Hall conductivity. One interesting direc- 
tion for future research would be to find ways to separate 
the two contributions to the intrinsic AHE in experiment, 
for example by measuring frequency or temperature de- 
pendence of the anomalous Hall conductivity. 
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